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Abstract 

We consider the fluctuation of linear eigenvalue statistics of random band n x n matrices whose 
entries have the form Mij = (\i — j\)wij with i.i.d. Wij possessing the (4 + e)th moment, 

where the function u has a finite support [—C*, C*], so that M has only 2(7*6+! nonzero diagonals. 
The parameter 6 (called the bandwidth) is assumed to grow with n in a way that b/n —>■ 0. Without 
any additional assumptions on the growth of 6 we prove CLT for linear eigenvalue statistics for 
a rather wide class of test functions. Thus we improve and generalize the results of the previous 
papers [8] and [11], where CLT was proven under the assumption n » b » . Moreover, 

we develop a method which allows to prove automatically the CLT for linear eigenvalue statistics 
of the smooth test functions for almost all classical models of random matrix theory: deformed 
Wigner and sample covariance matrices, sparse matrices, diluted random matrices, matrices with 
heavy tales etc. 


1 Introduction and main results 

Consider an ensemble of random symmetric n x n matrices with entries of the form 

M^J = Uij = u{\i - j\/b) (1.1) 

where are i.i.d. (up to the symmetry wij = Wji) random variables, satisfying the moment 

conditions 

E{wij} = 0, E{\wij\'^} = 1, E{\w^jf} = 3 + K4, E{\wij\‘^~^'^} < C < oo, (1.2) 

and u{x) is a piece-wise continuous (with a finite number of jumps) continuous at x = 0 function with 
a compact support, satisfying the conditions 

u{x) = u{—x), 0<u{x)<C, j u{x)dx = l, suppu C [—(7*,(7*]. (1.3) 

It is easy to see that the entries of M are nonzero only inside the band \i — j\ < C*b. Hence for fixed 
b we have a matrix with a finite numbers of diagonals, while if 6 ^ n, we obtain some kind of the 
Wigner matrix, with all of the entries having the variances of the same order (see m)- The model is 
now widely discussed in mathematical literatures, since by non rigorous conjecture of [5] it is expected 
that the behavior of local eigenvalue statistics demonstrates a kind of phase transition: for 6 << 
the statistics is of Poisson type and for 6 >> it is of the same type as for Wigner matrices. Till 
now this result is not proven rigorously, but the problem is one of the most challenging in the random 
matrix theory (see, e.g. m, la, a, m and references therein). 

It was proved many years ago (see [10]) that in the limit 

b oo, b/n —>-0, as n —i oo, (1-4) 

the normalized eigenvalue counting measure converges weakly to the Wigner semicircle low, which has 
the density 

PscW = ^V^4- A21[_2,2]- (1.5) 
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This means that if we denote the eigenvalues of A4, choose any bounded integrable test 

function (p, and consider the linear eigenvalue statistics of the form 

n 

NnW\- E{Nn[p>]}, (1.6) 


then in the limit m we have 

E{n~^J\fn[(p\} J ip{X)Pec WdX, Var{A/'„[v3]} 0. 

In particular, for </5(A) = (A — z)~^ 

n~^Afn[(p] = n“^Tr (M — z)~^ —>■ g{z), 

9{z) = i(-2: +(1.7) 

The next natural question is the behavior of the fluctuations M°['~p\ in the same limit, in particular, 
the behavior of its variance. This question was solved partially in the paper [5], where the main term 
of the covariance of the traces of two resolvents was found in the case of Gaussian Wij and under the 
additional restriction 6 = u®, 1/3<0<1. The next step was done in the papers [TT] and [5], where 
the Central Limit Theorem (CLT) for the random variable ^Jb/nN°[^p\ was proved for sufficiently 
smooth test functions, but again under the technical condition n » h» . 

The main result of the present paper is the proof of CLT for the linear eigenvalue statistics (II.6p of 
the band matrices under the limiting transition dUD without any additional restriction on the growth 
of h. 

We consider the test functions from the space Hs, possessing the norm 

ll<Plls = y(l + 2|A:|)^®|^(fc)pdA:, s > 2, ip{k) = J E^^ip{x)dx. (1.8) 


Theorem 1 Consider the model of band matrices and any test function possessing the 

norm with s > 2. Then the sequence of random variables y^6/nA/'°[v5] with of 

converges in distribution in the limit o to the normal random variable with zero mean and the 
variance 


V[lp\ =- 


dxdyip{2 cos x)ip{2 cos y) 


{u, u)k 4 


dxdy 
2 


log 


1 - u{k)e^^^+y'^ 


1 - u{k)ed^-y) 


dk 


(1.9) 


pTZ \ ^ T/fO) / f'^ 

J ip(2cosx)cos2xdxj + ^ ^ i^J (^(2 cos a;) cosxdx 


where (u, u) = f u^(x)dx and u(k) is the Fourier transform of the function u defined as in U.8\) 


To prove CLT for the band matrices, we use the CLT for martingale (see [U Theorem 35.12]). 


Theorem 2 Let Xn^k = E^k{Y — EkY} be a martingale differences array with respect to independent 
random vectors Vi,..., V„, Sn = ^k, = ^^=1 ^i^k} ~ *^(1)- Assume that 

n 

(1) (2) Var{^X2}<£„ (1.10) 

k^l 

Then 

\E{E^s^} - e-‘^-/2| < C\t){e]/^ + f^^). (1.11) 

Remark: Here we have replaced a more general condition ^ E{X^l^Xk\>s} 0 used in by condition 
(1) which is more easy to check for the random matrix models. 
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The idea to use Theorem [5] for the proof of CLT in the random matrix theory is not new. Since 
the paper of [T] it was used many times (see, e.g., [7], [E] and m), but the method of the proof 
of CLT used in the present paper allows to prove CLT by the same way for all classical models of 
random matrix theory: deformed Wigner and sample covariance matrices, sparse and diluted random 
matrices etc. It becomes even simpler than that for band matrices, since the proof of condition (2) 
becomes simpler. 

The paper is organized as follows. In Section 2.1 we give the sketch of the proof of CLT, introduce 
truncated band matrix and explain how one can extend CLT from some special class of the test 
functions to all functions of Hs- In Section 2.2 we check conditions (ll.IOl) and in Section 2.3 prove 
Lemma [T] (given in Section 2.1) about the difference of linear eigenvalue statistics of initial and 
truncated matrices. In Section 3 we compute the variance (HH). And in Section 4 the proofs of some 
auxiliary results (partially known before) are given in order to make the proof of Theorem [T] more self 
consistent. 


2 Proof of CLT 

2.1 Strategy of the proof 

We start from the proof of CLT for the truncated and ’’periodically continued” model: 

Mij = Uij = u{\i - j\n/b) (2.1) 

1 Wy, II* -j| -n\< C*b 

Here and below 

\i-j\n ■■= ma.x{\i - j\,\\i - j\ - n\}, ( 2 . 2 ) 

and {uJij}\\i-j\-n\<c*b are independent (up to the symmetry conditions) and independent from A4 
copies of Wi 2 - Thus we not only truncated the entries of A4, but also add entries in upper right and 
lower left parts of it, in order to obtain the periodic distribution, i.e., invariant with respect to the 
shift i ^\i + 1|„. 

Then the standard argument gives us that for \i — j\n < C*b 

E{w,j} = 0, A{|zcyf} = 1 + 0(6 -i-"/2), (2.3) 

} = 3 + «4 + 0(6-"/2), E{\LO,^f} < Cb^-^/^. 

Moreover, it is easy to see that 

n-^E{Tr (M - Mf} < Cb-^^'^. ( 2 . 4 ) 

Then, using Theorem[2 we prove CLT for vin := {b/ny^‘^Af°[ipn, M] with the test functions of the 
form 

‘Pr, = ‘P*'Pr„ ( 2 . 5 ) 

where * means a convolution, is a Poisson kernel 

( 2 - 6 ) 

and ip GTis Cl Li(R). It is easy to see that then 


Afn[Pr„M] = TT ^ 


(^(A)S7„(A + iri)dX. 


(2.7) 


Then we shall prove the lemma 
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Lemma 1 Set S(z) = (A4 — z) %(z) := Tr G{z). Then for any z :^z > t] 

^Var|7„(2;) - 7„(2:)| < (2.8) 

The lemma implies that for any Lp Li(IR) if we set v^n '■= then 


Var{i/2„ - vin} =^Var|A/'„[(/5^,X] - Mn[pr,, 


— J J dXidX2p{Xi)p{X2) 

X Cov{37„(Ai + iri) - 37„(Ai + i-q), 57„(A2 + iq) - Xs%{X 2 + iq)} 

< Cb-^ j J dXidX2\p{XiMX2)\ < 

Hence, for any fixed x S R 


Thus, CLT for vin and Lemma [T] imply CLT for V 2 m if the test function has the form (12.5|) . 

To extend CLT to the test functions from TLs, we use a proposition (see [Ill Proposition 3.2.9]). 

n 

Proposition 1 Let be a triangular array of random variables, A/"n[</5] = be its 

1=1 

linear statistics, corresponding to a test function (^ : R —>■ R, and {(i„} is some sequence of positive 
numbers. Assume that 

(a) there exists a vector space C endowed with a norm ||...|| and such that Vn is defined on C and 
admits the bound 

dnYe.r{Mn[p]} < CUvslP, Vp € C] (2.9) 

(b) there exists a dense linear manifold Ci G C such that CLT is valid for J\fn[p], p G Ci, i.e., there 
exists a continuous quadratic functional V : Ci ^ R+ such that we have uniformly in x, varying on 
any compact interval 

lim Zn[xp\ = e~^ "G[ip]/ 2 ^ where Zn[xp\ := (2.10) 


Then V admits a continuous extension to C and CLT is valid for all Mn[p\, p G C. 

The proposition allows to extend CLT from any dense subset of TLs for which we are able to prove 
CLT on the whole TLs, if we can check (1^ . This can be done by using the another proposition 
(proven in [16] and also [17]) and Lemma |2j 

Proposition 2 For any s > 0 and any A4 


poo poo 

Var{7V„[75, Ad]} < C's||(/j||g / Va.r{T t g{x + iy)}d: 

Jo J — oo 

Lemma 2 If the conditions E2P and are satisfied, then for any 0 < y < 1 

— f dxYailTTgLx + iy)} < Cy~'^\ogy~^ 
n J 


( 2 . 11 ) 


( 2 . 12 ) 


The proof of the lemma is given in Section 4. 

Combining the proposition with (I2.12L we prove (12.91) . 
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2.2 Checking of conditions (ll.lOp 

To apply Theorem [21 we denote Ep the averaging with respect to the variable {wp,j}j>p, E^p = 
El... Ep_i and consider 


XpYpr,]=T^-^{b/nf/^ j ip{\)<^Xp[\ + iri\d\ 


(2.13) 


Then, according to Theorem|2l we have to check condition (l)-(2) of p.lOll for {Xp[^pp\}. It is evident, 
that conditions ( 1 ) follow from the bounds 

Ep{\Xp[z\\^}<Cb-\ \Xp[z]\<C, (2.14) 

valid uniformly in |5Jz| > 77 . And since 

Var{ E - E Ep{Xl}] = E I^P - Ep{Xl}\^) < E E{X^p} < 

p p 

condition (2) of (I1.10|) follows from the uniform in |3zi|, |322| > ?? bound 

Var{S(zi,Z2)} <en, (2.15) 

S(zi,Z2) := ^ E Ep{Xp[zi]Xp[z2]}. 


Let us prove (12.141) and (12.1511 . 

Denote M^Ei the ( n — 1) X (n — 1) matrix which is obtained from M by removing the pth line and 
column. Set also 

G(P) = (M(p) - z)-i, := (upi,... ,Up„) e v,, := (2.16) 

Use the identities 

Gpp = -Ap\ G,i=Gif-Q[f, TtG-TyG^p^ = -^JogAp{z), (2.17) 

where 

Ap := z + b-E\^ + b-^(G(P^v(P\v^P^), (2.18) 

Qlf = b-^Ap^{G‘^P'>v^P'>UG<'P^v^P^)j. 

Since for the resolvent G{z) = (M — z)~^ of any symmetric or hermitian matrix M and any vector m 

’A(G{z)m,'m) = 3z(G(z)m, G(z)m), (2-19) 


we have for |3z| > rj 


and 


\Ap{z)\ > |3Ap(z)| = |3z|(l + 6 -i(G(p)n(^’\G(^’)'i;(p))) > 77 , 
|Ap| > \?sAp\ > 77 , where Ap := Ep{Ap}, 

IIQ^^^II < \Ap\-^\b-^{G^P'>v^P\G^P'>v^P'>)\ < 77 


A'piz) |i + 5-i((G(p))Vrf, 


Ar, 


< 






|Ap| < 2,- 


which implies the second inequality of (12.141) . 


( 2 . 20 ) 

( 2 . 21 ) 
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The last relation of (I2.17P yields 


Ep{Xp{zi)Xp{z2)} =■■ Dp{zi,Z2) 

Dp{zi,Z2) ■■= Ep^E<p{{logAp{zi))p}E<p{{logAp{z2))°p}Y 

Here and below for any random variable ^ we denote Cp = ^ ~ ^piO- 

Since Dp{zi,Z2) is an analytic function on zi,Z2 '■ I O'- 22 ! > ??/2, in order to prove the first 

bound of (12.141 it suffices to prove that uniformly in > ??/2 

Ep{\E^p{i\ogAp{z)rp}\^} < r^Ep{\E^p{Al{z)}\^} < Ch-\ 


Evidently 

E<p{A°p(z)} = b-^/^Vpp + b-^ G[f{z)vpiVpj + b-^ ^ (z)(Upi - Upi). 

i>P 

Hence, averaging with respect to Ep and using (12.31) . we obtain the first bound of (I2.14p . Similarly 
one can get the relation which we need below 

Ep{\E^p{Al{z')}\'^} < Cb-^--/\ (2.22) 

We are left to check (12.151) . Writing Ap = Ap + Ap, expanding logHp around Ap, and using (12.221) . 
we obtain 

S(zi, Z 2 ) ^ Ep{Xp{zi)Xp{z2)} = Z2) (2.23) 

E{zi,Z 2) =:^^Dp{zi,Z2) = ^^{Ap{zi)Ap{z2))~^Tp{zi,Z2) 

+ ^ E {0{Ep{\A;{z^)\^}) + OiEp{\A;{z 2 )f}) 

=- J2{Ap{z,)Ap{z2))-%{z,,z2) + 0{b-^/% 

where 

Tp(zi,Z2) ■.=Ep{E^p{Al{z,)}E^p{Al{z2)}} (2.24) 

=26-2 ^ Up.Up,Ep{E^p{G‘'^\z^)}E^p{G‘f^{z,)}} 

i,j>P 

+ AC 46 - 2 ^u 2 ,Ep{E<p{G(f(zi)}£;<p{G(f)(z 2 )}}+ 6 -'^Xpp. 

Lemma 3 Given ry > 0 there exists S(r]) > 0 such that uniformly in z : |52:| >77 

Var{Gg)} < 6-^ E{\G^^ - E{G,,}\} < Clb-\ (2.25) 

E{\G^^(z)-g{z)\}<C^b-\ 

where g{z) is defined by (Z3- 

The proof of the lemma is given in Section 4. 

Remark 1 Below we will often use a simple observation. If for some random variables |i?fc| < Ck, 
Sfc C'fe < C, and fk ■ E{\fk — fk\} ^ C'l^ where ff. are some constants, then we have with the 
same G and Gq of i2.25\) 


Y,Rkfk = J2Rkrk+r, E{\r\}<GC,b-^. 


(2.26) 
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(2.27) 


In particular, since in view of (12.241) \Tp{zi,Z 2 )\ < Cb we have 

f^{zi,Z2) = ^'^{Ap{zi)Ap{z2))~^Tp{zi, Z2) + 0 ( 1 ) 

2 

= -^'^ 9 {zi) 9 iz 2 )Tp{zi,Z 2 ) + K4{g{zi)giz2)f + u{ 0 )g{zi)g{z 2 ) + o{l), 

where Tp{zi,Z 2 ) is the first sum in the r.h.s. of (12.2411 . The constant term here does not contribute 
into the variance of S( 2 i, 2 ^ 2 ), so it is not important in the proof of (12.151) . 

Let us denote the matrix M whose entries Wij with min{j, j} < p are replaced by Wij which 

are independent from all {wki}^i^i and have the same distribution as Wij. Let also be the 

matrix M^^p^ without qth line and column. We denote also E^p the averaging with respect to all Wij 
and Wij with min{i, j} < p. Set 

G«P’9) = - z)-\ - z)-\ (2.28) 

Then evidently 

T;{zuZ2) = (^2)}}zi,pUfcp 

jk 

= Ep{E<p{TTG^<P’P\zi)I^P^G^P\z2)I^P^}}, 
where we denote by I^p^ the diagonal matrix with the entries 

= Sj,kUkplk>p- (2.29) 

Moreover, if we replace G^p'^ in (12.241) by G and set 

Tp{zi,Z 2 )= ^ UpiUpjEp{E<^p{Gij{zi)}E<^p{Gij{zi)}} (2.30) 

i,i>P 

=Ep{E^p{i:T&<p'> (21) G(22)/(p) }}, 

then in view of (12.171) and (12.211) 

|T;(2 i,22 ) -T;(2 i,22 )| <|Gp{G<p{TrQ(P)(zi)/(J’)G(22)/(P)}}| (2.31) 

+ \Ep{E<p{TY&<P^P\zi)l‘^P'>Q^Ph^P^}}\ < G, 

where we have used that since Q^p^ is a rank one matrix with a bounded norm, we have for any 
bounded matrix B 

TrQ^P'>B<\\B\\\\Q\\. 

Thus we need to study the variance of 

(2.32) 

V 

To prove ()2.15p . it suffices to show that 

Var{Ei} = ^ EilElAi^.yX} ^ 0. 

r 

The last relation is a corollary of of the bounds, which we are going to prove 

n2G{|(Ei)°|2}<G, r = l,...,n. (2.33) 

By (EJll), 

= (2.34) 
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Notice also that (r"(zi, 2 : 2 ))° = 0 for p > r + 1, hence the sum in (12.351) is over p <r. 

Then (12.171) yields 

(r;(zi,22))° =(^<p{TrG«P)(zi)/(P)G(22)/^^)} - ^<p{TrG«P’’')(2i)/(P)G«(22)/(^H) 
= {e<p {{Arb)-\G^^^ {Z2)I^P^ G^<P’A (zi)/(p)g(’’) iz2)v‘^^^ , ) ly 

+ Sim + (g<p {(A, 6) -2 (G(") (22) G(") (zi )2 } y 

=:{F^p\zi,Z2))l + {F^p\z2,zi))l + {F^p\zi,Z2))l, 


where ” +sim” means the adding of the term which can be obtained from the previous one by replacing 
Z 2 and Zi- Since < -G{|CP} for any random variable p.34l) yields 


i?{|(Sl);P} < Gg{|6-1^ (.F«(zi, 22) +Fi;Hzi,^2))f } 

p<r 

< Gg{ |6-2 T;<p{(/(p)g(’')(zi)u(’'\ GM(zi)u(’')) (l + b-\v^A } j^} + sim 

p<r 

=: GT;{|&-2^ (fM(2i) +FM(22))f }. 


(2.35) 


p<r 


To sum in the r.h.s of (12.351) with respect to p we would like to use the property 


Y < Cbl, (2.36) 

p^l 

but since p appears not only in I^p\ we need to remove p from the other places first. Write 

E{n^\{j:YS} <Cb-^ E 

p<q<r 

r q 

<G6-4 EEM E<q { (/(P) G^AyG) ^ Q(r)^(r) ) (1 + 6-1 (^(O ^ ^ W ) )jp(r) | 

9=1 P=1 
r 

<G 6-2 E ^<9 { ) (1 + } ^ 9 ’’^} 

9=1 

r—C+6 

9=1 

+ CE{{l + b-\v^A^y(r)))4-^ < c'T;{(1 + 6-1(z;("),u(")))4}. 

Here in the first line we use (I2.35P , in the second line we use first that for p < q the averaging E<p can 
be replaced by E<q, and then use (I2.36|) for summation over p < r. The third line follows from the 
second one in view of the bound ||G^'’^|| < G. Next we split the sum over q into two parts: one over 
q < r — G*b and another over r — G*b < q < r, and observed that for the q in the hrst part {v^'^\yir)'j 
is a constant with respect to the averaging E^q, hence 

E{E^q{b-\v^^\v‘^^^){l + b-\v^^\v^^^))}E^A^ 

= G{((GM)*/(«)g(’’)uW,u(’'))5-1(u(’’\u(’'))(1 + 6-1(uW,uM))^}. 

Then we can take the sum over q < r — G*b, using again the bound (12.361) . and finish to estimate the 
sum using the bound ||G*^’'^|| < G. As for the terms with r — C*b < q < r, they are estimated just 
using the boundedness of ||G^’'^|| and Thus we have proved (12.331) . 

□ 



















2.3 Proof of Lemma [T] 

Set 

g(P^ := (M^P^ - z)-\ Ap:= z + b-^{g^P^d(P\d^P^), AAp := Ap - Ap. 
The same argument as in the previous section implies that it suffices to check that 

- ^ E{\AAp - Ep{AAp}\'^} ^ 0. 


Since we know that (see (12.241) 1 
b 
n 


^ E{\AAp-Ep{AAp}\^}<- Y. 2(^E{\A;n + E{\A;n)<—, 

\p\n<C*b ‘ bU<c*& 


we conclude that it suffices to prove that 

- V E{\AAp - Ep{AAp}\^} ^ 0. 

n ^' 

|pU>c*b 

Let us write 

AAp =b-^/^Avpp + b-^{g^P^ Av^P'> , Av^P '>) + 2b-\g^P^ Av^p^ ,v^p^) 
+ b-\{g^P'> - G^P^)v^P\v^P^) =: Jop + Jip + 2J2p + hp. 

Averaging with respect to and v^p^ we get similarly to (12.241) for \p\n > cb 

E{\Jip-Ep{.hp}\'^} = b-^YE{\Q\f\'^('^pi - - ^Pjf} 


U-l-e/2 


|2}/W/W +6-2-e/2^^{|g(P)|2|^W < Cb-^- 

i^j i 

Similarly 

E{\J2p - Ep{J2pW} < E{\Jop - Ep{JopW} < 

In addition, again similarly to (12.241) we have 

^^{|J3p - ^ip{^3p}p} <cb-'^E{TY&\g'^p'^ - G^p^)i^p\g^p> - G^P)*)}. 

Now by the same way as in (I2.30I) - (I2.31I) we can replace here g^P'^ by g and G^p'^ by G with 

0 {b-^): 

E{\J3p - Ep{J3p}\^} <2b-^E{TTl^PHg - G)I^P\g* - G*)} + 0(6-2). 

The resolvent identity implies 

g-G = G{M^P'> - M)g = -GAMg. 

Hence, the last term in the r.h.s. of (I2.42p can be estimated as 

b-^EiTiI^P^g - G)&\g* - G*)} = b-"^E{T t I^P'^GAMg&^g* AMG*)} 

< Gh-‘^E{Ar&'iG{AMfG*)}. 

Hence, using (12.361) and (12.41) . we obtain 
b 
n 


Y E{\J3p - Ep{J3p}\^} < Gn-i6-iA{TrG(AM)2G*} 

C*b<p<in—C*b 

< Gn-^b-^E{TT{AMf} < Gb-^-^/^. 


Combining (12.441) with (12.391) - (12.421) . we get (12.381) . 

□ 


(2.37) 

(2.38) 

(2.39) 

(2.40) 

(2.41) 

(2.42) 
L error 

(2.43) 

(2.44) 
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3 Variance 


In view of (12.321) to find Si, it suffices to find the main order of b ^£'{T"(zi, 2 : 2 )} defined in (12.301) . 
For this aim it suffices to compute for any i the main order of 

ti = ^^Upj E^p{Gij (^Zi^Gij(^Z2^\- 
3>P 

Consider 

Si ■■^'^UpjEcp^Gij{zi)'^h~^/'^VikGkj{z2)^ (3.1) 

3>P k 

— ^ ' upjE^p^Gij{zi) ^^^ ^ Vik Z2Sik C/cj (^ 2 ) ^ 

3>P k 

= X! UpjSijE{Gu{zi)} + Z 2 ti = Upig{zi) + Z 2 ti + 0{b~G'^)^ 

3 

where we used Lemma |3] for the last equality. 

The idea is to compute the l.h.s. above in a way which gives us an equation with respect to {ti}i>p. 
It is possible by using the formula (see e.g. El) valid for any random variable ^ which has zero mean 
and possesses m + 2 moments, and any function F, possessing m + 1 bounded derivatives 

E{mO} = E + r^, |r^| <CF;{|er+nmax|i^(™+i)|. (3.2) 

Applying this formula for ^ m = 4, and Fijk = Gij{zi)Gik{z 2 ), we get 

Si = ^ ', ^pjE<Cp^Gji(^Zi)Gij(^Z2) ^ ^ b UikGkk(^Z2)^ ( 1 ^- 1 ^) 

3>P k 

- ^ 5 ^u^kUpjE^p{Gu{zi)Gjk{zi)Gkj{z2) ^ Gjk{z2)} 

3>P k 

+ Ri + R2 + R3 + Ri- 

Here we used the differentiation formula for the resolvent of any symmetric matrix M 

-^Gsi{z) = -Gsk{z)Gu{z) - Gs^{z)Gkl{z) (3.4) 

dMik 

Two sums written in the r.h.s of (13.31) collect the terms, corresponding to s = 1 in the r.h.s. of (13.41) . 
The remainder i?i collects the terms, corresponding to s = 2 in the r.h.s. of p.4p . The remainders 
i ?2 and i ?3 collect the terms, corresponding to s = 3 and s = 4 respectively. And the remainder R 4 
appears because of the remainder in (13.21) . Let us analyze the order of each of these terms. By (13.41) 

Rl — b ^ ( UjpUikE^p\^Gij(^Zi')Gik{z2)Gkj{z2^\ 

3>P,k 

b ^ ( y'jpy'ikE<^p\^Gik{zi)Gij(^Zi)Gkj(^Z2^\ 
j,k>p 

= - b-^E<p{{GI^P'>GI^^’P^G)^^} - = 0(&-i). 

where = 6 ikUiklk>p- 

To estimate R2, observe that by (13.41) after two differentiation we obtain the sum of terms of the 
type Gi-^i^Gi^i^Gi^igGijjg, where G can be G or G and the set of indexes Zi, I 2 ■ ■ ■ I 7 , h contains 3 times 
f, 3 times k, and 2 times j, but Gjj can not appear. Thus, each term contains either GjkGji or 
GjkGjk, or GjiGji- Any of this combinations after summation with respect to j gives us 0(1). Hence, 
after summation with respect to k we obtain 0{b). But the factor which appears because of the third 
cumulant is b~^G^ hence i ?2 = 0(6“^/^). By the same argument R^ = 0{b~^). 
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Finally, to estimate i? 4 , observe that we have two summations with respect to p < j < p + C^:b and 
i — C^b <k <i + C^b, and the factor which appears because of b~^E{\vik\^} is bounded by 
At the last step of transformations of (13.31) we write 


Gkk{z2) = g{z2) + iGkk{z2) - g{z2)), Gii{zi) = g{zi) + {Gu^zi) - g{zi)) 
and use the bound (12.251) . Then we obtain 

Si = -g{z2)ti - g(zi) ^ Ul^hk + ri, u < Cb~G'^^ 

k 

where 

Uik — b Uik-, G — b Uik^i>p^k>p' (^•^) 

Combining m and (1331) with above estimates for the reminders and using that by dni we have 
{z 2 + giz 2 )) = —g~^iz 2 ), we obtain the system of equations 

((C - + r', r' < C{b-G^ + b-G^), 

with C = {g{zi)g(z 2 ))~^, = li>pUpi. 

Since \g{zi)g{z 2 )\ < 1 and 

||C/^^)|| < max^ \Uki\ < 1 + o(l), 

i 

the operator (C — can be defined by the Neumann series 

OO 

(C - ^ ^ 

m—0 

and it possesses the properties 


(3.6) 


-0-k\<G, i>p, \{U'^U<Gb-\l<i<n. (3.7) 

k 

Application of {( — to both parts of (13.61) and (13.71) imply 

t^=([U^p) + \r^\<Cib-G^ + b-G^), (3.8) 

^b-%'izi,z2) = 6-1 ((C - + o(l), 

^A{Si} = A^6-i((^_c/(rt)-i>),y(p)) +o(l), 

^ p 

where Si was dehned in (|2.32p . 

Proposition 3 Let the matrices U and he defined by SS.5\} . where {uij} satisfy conditions (LM>, 
the vectors jg defined by hS.bft . and |C| > 1. Then 


1 h 

-^y"b-\{C-U^Gyiuip)^uip)) = --fTrlog(l-C-if/) FC-iTrc) +0(6-1). 

p^i 

Proof. Denoting by S'i( 2 :) the l.h.s. of (13.9p and by S 2 {z) the main term in the r.h.s. of ()3.9I) . we have 

, OO 

S2iz) =- V m-iC-™ V 

n 

m—2 

,7100 

X! X! 1” ™ X! ^* 1*2 • • • 

p—1 m—2 min{ii,...,i,n}=P 

, 71 OO 

=-EEc-™ E +0(6-1) 

p—1 m—2 i 2 ,...,im-i>P 

^ 71 OO 

=^EE ^ ^ ^ 0(6-1). 

p—1 m—2 


11 






The term 0{b~^) appears in the third line above as a sum of the terms, which have at least two p 
among {ii ,..., im\- But the contribution of these terms for fixed m in view of (13.7p can be estimated 
as 

m—1 

fe=i 

After summation with respect to m and multiplication by b we obtain 0{b~^). 

□ 


Now observe that the r.h.s. of (13.911 has a limit, as n, 6 —>■ oo like in (11.41) . 

CXD 

S2{Z) = 

m—2 


-i— / u{xi - X2) ■■■ u{Xm-l - Xm){u{Xm) “ u{xi))dx + 

J 


= —J\og{l-( ^u{k))dk-( ^u( 0 ) + o(l). 


where u is the Fourier transform of the function u defined as in HH). Hence, the proposition and the 
last line of (13. 8|) yield 


F;{Si} = -2 (y log (1 - C^u{k))dk + C“^u(0)) + o(l). 

Thus by (12.231) and (12.271) we obtain 

^Cov{7(zi),7(z2)} = E{T,(zi,Z2)} 

= ( ~ 2 J log {1 - g{zi)g{z 2 )uik))dk - g(zi)g{z 2 )u{ 0 ) + K4g'^{zi)g'^{z2)^ + o(l) 

=: C{zi,Z 2 ) + 0 ( 1 ), 

where we used also that by (13.61) = g{zi)g{z 2 )- According to the definition (12.71) and the above 

relation 

b f 

-Vdx{Nn[Pr^} / d\idX2ip{\l)ip{\l)Cr,{Xl,\2), 

where 

^7,(^1, A 2 ) =^(C(Ai + Z77, A 2 - ig) +C(Ai - zt?, A 2 + ig) 

— C{Xi + ig, X 2 + ig) — C(Ai — ig, X 2 — ig)). 

Now by Proposition [ 1 ] for any p possessing the norm (11.81) we have 

b f 

-Var{A/'n[v3]} -t lim / dAidA 2 (/?(Ai)(p(Ai)C'^(Ai, A 2 ). 
n J 

Let us make the change of variables Ai = 2cosa::i, A 2 = 2 cos 0 : 2 . Then, using that (see dm) 

lim g{Xadszig) =, a = 1,2, 

17-*-+0 

we obtain dm by a simple calculus. 


4 Auxiliary results 

Proof of Lemma 

The first identity of (j2.17ll yields that it suffices to estimate E{\ApA~^ — Ei{A'pA~^}\'^}. Note 
that for any a independent of {wii} we have 

Ep{\Q^}<Ep{\i-a\^}. 
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Hence it suffices to estimate 


■^p Ep{Ap} 


4'o 

/ip 

A'° A' 

-^p -^p 

< 

A'° 

/ip 


A°p 



Ep{A} 

E'i{Ap} A.p 


Ep{Ap} 


yEp^Ap} 


Here and below z = x + iy, y > 0. Let us use also the relation (I2.19|) which yields, in particular, that 
\^p/Ap\ < y~^ . Using (12.181) , we get 


Ep 


II 


Ep{Ap] 


< 


Cb-^Tv G^P'flpG'^Ph 

\Ep{Ap}\^ ’ 


Similarly 


Thus 


Er. 


A'° 

/ip 


Ep{Ap} 


< 


Cb-^Tr {G^P'>flp{G^P> f Cb-^Tr G^P'>IpG^P'> 


\Ep{Ap}\^ 


< 


\Ep{Ap}\^ 


^ _iv-TrG(P)/pG(P> 
-E{\hn{z)) I }<Gn ^ by^\Ep{Ap}\'^ ' 


Notice that the Holder inequality implies for any (5 > 0 

J\b-^ H gI^\x + iy)\^^^dx < Gb^^ J \G^/^\x + iy)\^+^dx 


\j-p\<bC 




|(cC-Afe)2+y2|(l+5)/2 


\j-p\<bC" 
dx < GS-^y-^. 


Hence, denoting Cp = {x ■.\^ UpjG^^^ (x + iy)\ > 1}, we obtain 

J ^Cpdx < Gmm{(5“^y“^} < Glogy” 

Then, using once more that by (12.191) each summand in the r.h.s. of 
/ -E{\{-fn{zx + iy))°\‘^}dx 


(4.1) 


(4.2) 


is bounded by y we get 


< Gn ^ ( [ iy^b) G{x + iy)IpG^P'’* {x + iy)dx + Gy ^ f dx^ 

^dR\([-l,l]u£p) Jl-l,l]UCp ' 


<Gy^yGy "‘logy ^<G'y '‘logy ‘. 


□ 

Proof of Lemma 0 It follows from (12.171) that 

E{\Gpp - E{Gpp}\^} < \Qz\-^E{\Ap - E{Ap}\^} (4.3) 

< 2\<Az\-\E{\Ap - E{Ap}f} + Var{H;}) 

< 2\<Az\-%-^'^Up^ViiT{Gif} + Gb-\ 

But since 

E{\Gu-GI^^\} < \^z\-%-^E{\{G^p'^vG)),\^} < G\^z\-%-\ 

we have 

Var{G,(f } = Var{G,J + 0(&-‘) = Var{Gpp} + 0(6-‘). (4.4) 

Here the last equality is due to the invariance of the distribution of M with respect to the ’’shift” 
i —>■(* + 1) mod (n). Hence for any 2 : |5>z| > 2 we obtain from (14.31) 

Var{Gpp} < 2|5z|“^Var{Gpp} + G6“‘Var{Gpp} < 2Gb~^. (4.5) 
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Let us fix any z = x + ir] with 0 < rj < 2 and consider the function 


0(C) = log(co6^/^|Cov{C?pp(C),Gpp(z)}|) 

in the half-circle 17 = {SC < 2} n {|C — x — 2i\ <\2 — r]/2\]. It is a harmonic function, and in view of 
(I4.51) for SC = 2 we can choose cq sufficiently small to have 

co6i/2|Cov{Gpp(C),Gpp(z)}| < co6i/2Vari/2{Gpp(C)}Vari/2{Gpp(z)} < 1 

^0(C)<O, C e 71 := 917 n (SC = 2}. 

Moreover, in view of the trivial bound |Gpp(C)| < |SC|“^, we have 

0(C) < log 6 ^/^ -l-logcory"^, C e 72 := 917 n {|C - a: - 2i| = |2 - ry/2|}. 

Hence, by the theorem on two constants (see [S], p. 296), we have 

0(C) < (log&^/^+ logco?7"^)w(C), (4.6) 


where the harmonic function 


w(C) := 



2 - 77/2 - (C - a: - 2 i) 
2-r7/2-h(C-a;-2i)’ 


satisfy the conditions 

w(C) = 0, C e 71 , a;(C) = 1, C e 72 - 

Since u]{z) = 1 — 25 with some 5{r]) > 0, (14.61) implies the first line of (|2.25|) : 

co&^/^Var{Gpp( 0 )} < ^ Var{Gpp(z)} < G&"0 


Using (I2.17L (j4.3l) . and (14.6L we get similarly to (l4.4L 

E{Gpp{z)} = -{z + E{Gpp{z)})-^ + Oib-^) 
^E{Gpp{z)} = g{z) + 0{h-^). 


Thus, we have proved the second line of (I2.25p . 

□ 
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